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ABSTRACT 


The  theory  group  at  AFRL/PRSP  is  engaged  in  an  ongoing  effort  to  develop,  implement,  and  test 
the  spectral  theory  of  Schrodinger  eigenstates  for  efficient  construction  of  accurate  many-body 
potential  energy  surfaces  for  HEDM  systems  such  as  atoms  embedded  in  solid  hydrogen. 

Because  good  experimental  data  is  available  for  AIAtn  clusters,  they  make  an  excellent  test  bed 
for  the  spectral  theory  and  its  computational  implementation.  We  have  used  the  spectral  theory 
method  to  calculate  the  ground  and  low-lying  excited  states  of  AIAtn  clusters,  subsequently 
employing  these  in  a  sequence  of  quantum  and  classical  Monte  Carlo  and  classical  molecular 
dynamics  simulations  to  predict  their  structures,  photoionization-induced  fragmentation  pathways, 
and  optical  absorption  spectra.  These  simulations  address  the  key  issue  of  the  degree  of  cluster 
fragmentation  encountered  following  photoionization  of  the  neutral  clusters.  The  measured  data 
refer  to  ion  mass  spectrometric  signals  of  [A1Atn-m]+  fragments  detected  subsequent  to  one- 
photon  uv  excitation  followed  by  visible  one-photon  ionization  of  a  distribution  of  non  size- 
selected  parent  AIAtn  clusters.  Thus,  the  observed  spectrum  attributed  to  AIAtn-m  is^  actually  a 
superposition  of  spectra  of  all  parent  neutral  clusters  AIAtn  which  produce  [AIAtn-m]  as  one  of 
the  ionic  fragments  following  the  photoionization  step.  The  simulations  involve  three  steps.  (1) 
generation  of  the  ground  state  wavefiinction  vj /  of  the  parent  neutral  AIAtn  cluster  using  diffusion 
quantum  Monte  Carlo  methods,  (2)  vertical  excitation  of  an  ensemble  of  configurations  drawn 
from  the  ground  state  probability  density  |vj/|2  to  the  [AIAtn]  ionic  surface,  followed  by  classical 
molecular  dynamics  simulations  of  the  relaxation  and  fragmentation  processes  leading  to 
formation  of  a  distribution  of  fragment  ions  [A1Atn-m]+  +  M  Ar;  M=l,2,3,....,  and  (3)  classical 
Monte  Carlo  simulations  of  a  “hybrid”  absorption  spectrum  as  a  weighted  average  of  the 
predicted  spectra  of  the  neutral  clusters  AIAtn  contributing  to  the  formation  of  the  ionic 
[AlArN.M]+  fragment. 

f  RDL  Summer  Faculty  Research  Associate,  1998. 
tf  NRC  Senior  Research  Associate,  1998/1999. 


GOAL 


Stabilize  up  to  5%  atoms  in  solid  hydrogen. 

OBJECTIVES 

Develop,  implement,  and  apply  efficient  computational 
methods  for  reliable  prediction  of  properties  (e.g., 
structures,  stabilities,  electronic  and  vibrational  spectra)  of 
doped  cryogenic  solids. 

APPROACH 

Develop  efficient  computational  methods  for  reliable 
modelling  and  simulation  of  advanced  cryogenic 
propellants  such  as  boron  atoms  in  solid  hydrogen, 
including  (a)  construction  of  accurate  potential  energy 

surfaces  using  the  spectral  theory  of  Schrodinger 

eigenstates,  (b)  reliable  methods  for  proper  treatment  of 
quantum  effects  such  as  tunnelling  and  zero-point  energies 
(e.g.^diffusion  quantum  Monte  Carlo,  path  integral  Monte 
Carlo,  Centroid  Molecular  Dynamics),  and  (c) 
development  of  models  for  reliable  prediction  of 
laboratory  observables  (electronic  absorption  and 
emission,  vibrational  spectra^) 


OVERVIEW 


Although  AlArN  clusters  are  not  HEDM,  they  are 
nonetheless  useful  prototypes  for  actual  HEDM  systems 
such  as  B(H2)n  clusters  and  solids.  Because  good 
experimental  data  is  available  for  AIAtn  clusters 
(including  information  on  ground  and  excited  state 
potential  energy  surfaces),  they  make  an  excellent  test  bed 
for  our  new  spectral  theory  of  Schrodinger  eigenstates  and 
its  computational  implementation. 

EXPERIMENTS 

AlArN  clusters  are  generated  by  laser  vaporization  of 
aluminum  metal  and  supersonic  expansion  of  an  Ar  carrier 
gas.  The  clusters  are  then  ionized  by  1+1’  two-color 
multiphoton  ionization,  first  exciting  the  A1  atom  3p  ->  3d 
transition  with  the  doubled  UV  photon  and  subsequently 
ionizing  with  the  visible  fundamental  photon.  Following 
ioniziation,  the  clusters  are  detected  using  TOF  mass 
spectrometry. 
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SIMULATION  PROTOCOL 


1 .  DQMC:  Ground  state  wavefunction  of  neutral  AlArN 
parent  cluster  is  represented  by  collection  of  500 
“walkers”,  which  are  equilibrated  for  10-20  ps. 

2.  Each  “walker”  provides  an  initial  configuration  for  an 
MD  trajectory  on  the  ionic  potential  energy  surface. 
Vertical  excitation  from  ground  state  surface  to  ionic 
surface  is  followed  by  250  ps  trajectory. 

3.  The  size  of  the  cluster  at  the  end  of  each  MD  trajectory 
in  step  2  is  weighted  by  the  proper  descendant  weight  of 
the  originating  walker.  Results  are  averaged  over  all 
walkers  to  obtain  predicted  fragmentation  pattern. 

4.  Absorption  spectra  of  each  neutral  cluster  AlArN  which 
produces  a  fragment  of  size  AlArN_M  is  simulated  using 
AlArN  ground  and  excited  state  potentials  constructed  by 
spectral  theory  method,  used  in  conjunction  with  classical 
Monte  Carlo  simulations  for  ensemble  averages  of  AlArN 
cluster  geometries. 

5.  A  “hybrid”  spectrum  is  constructed  as  a  weighted 
average  of  the  spectra  of  the  parent  clusters,  using  the 
distribution  of  fragment  sizes  obtained  in  step  3 .  Uniform 
distribution  of  parent  neutral  clusters  is  assumed! 


SPECTRAL  THEORY  OUTLINE 


Diatomic  wavefunctions  written  as  superposition  of 
atomic  product  states  (direct  product  of  A1  and  Ar  atomic 
states): 

A1  Ar 

{3P,  4S,  3D,  4P}  ®  {ISO} 


®  =  {  3P+1.'S0,  3P0»iSo,3P.i»1S0,4S»1So,3D+2»1So,...,4P.i«1So} 
=  { 3P+1,  3P0,3P.1,4S,3D+2,...,4P.1} 

=  { x  2n,  a2s,  b22,  (3)  2n } 


d>  =  u,vF 

Hamiltonian  for  atomic  cluster  given  by: 

H  =  2  H(k)  +  H(SO) 

H(k)=  D‘(0k^k)  U‘(Rk)  E(Rk)  U(Rk)  D(ek,())k)  -E(#) 

H(so)  =  spin-orbit  coupling  operator  for  A1  atom. 

E  =  diagonal  matrix  of  diatomic  interaction 
energies  (X,A,B,C,...  state  energies) 


U  =  unitary  matrix  which  transforms  from  the 
diatomic  basis  to  the  atomic-product  basis. 

D  =  transformation  matrix  connecting  the  “rotated” 
and  laboratory  coordinate  systems. 

Atomic  state  mixing  parameters  Uij  obtained  from 
eigenvectors  of  the  diatomic  and  “spectral-product” 
transition  moment  matrices:  U  =  Ud  (U^)1,  where  Ud  and 
Usp  obey  the  following: 

A*  =  (Ud)1  m)  Ud 

and 


=  (Uspy  n»(R  -->  oo)  Usp 
If  Xsp  =  Xd,  then 

(udy  ii(r>  ud= (uspy  n«usp 

(Usp)  (Udy  n(R)  Ud  (Usp)‘  =  n* 

U'  n(R)  U  =  Ho 

This  is  “exact”  only  in  the  case  where  A,sp  =  Xd. 

U  is  related  to  diatomic  nonadiabatic  coupling  matrix 
elements.  If  Py(R)  =  <Ti  |  d/dR  |  Tj>,  then 

dU/dR  +  PU  =  0  =>  P  =  -  (dU/dR)  U‘ 

(Werner  and  Meyer,  JCP  74,  5802  (1981)) 
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CONCLUSIONS 


1.  A1  is  located  on  surface  of  clusters.  Anisotropy  of 
Al-Ar  ground  state  potential  energy  surface  favors  surface 
site  for  A1  atom. 

2.  Several  Ar  atoms  (2-13)  evaporate  from  the  neutral 
parent  AlArN  clusters  following  photoionization. 

3.  The  observed  and  simulated  spectra  assigned  to  the 
AlAri2  cluster  are  in  reasonably  good  agreement.  Less 
satisfactory  agreement  between  theory  and  experiment  is 
obtained  for  AlAr6.  This  discrepancy  may  be  due  to  the 
assumption  of  equal  populations  of  the  parent  neutral 
clusters  AlArN. 

4.  Spectral  theory  method  is  a  computationally  efficient 
method  for  calculation  of  reliable  ground  and  excited  state 
condensed  phase  potential  energy  surfaces.  Spectral 
theory  properly  accounts  for  avoided  crossings  of 
adiabatic  potentials. 

5.  Our  model  for  predicting  electronic  absorption  spectra 
of  atom-doped  rare  gas  solids  and  clusters  shows  overall 
satisfactory  agreement  with  experiment. 


FUTURE  DIRECTIONS 


1.  Simulate  relaxation  on  the  [AlArN]  intermediate 
excited  state  surface(s)  prior  to  excitation  to  the  ionic 
state. 

2.  Incorporate  non-additivity  into  ground  state  potential 
energy  surfaces  (both  Al-Ar  and  Ar-Ar  non-additive 
contributions^ 

3.  More  accurate  treatment  of  spin-orbit  effects 
(recalculate  diatomic  potential  curves,  including  SOCfjj 

4.  Combine  spectral  theory  potential  energy  surfaces  with 
quantum  simulation  method  (diffusion  MC,  CME^j 
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